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1. Introduction 

In recent months there has been a revival of interest in the noncritical c = 1 bosonic 
string |], 0, § (and its worldsheet supersymmetric counterparts 0, §,f§). New insight 
has been gained into this relatively simple string theory, using several developments 
that came after the previous matrix revolution: D-branes[0], M(atrix) Theory^, [|, 
AdS/CFTjll]] and the understanding of boundary states in Liouville theory [|TT|, [12], 



Most of the recent work has centred on matrix quantum mechanics, whose 
double-scaled limit is believed to represent the c = 1 string (at least perturbatively). 
The basic idea is that the matrix of this model is the matrix-valued tachyon on the 
worldline of a collection of N DO-branes. 

There are also claims that a nonperturbatively consistent version of matrix quan- 
tum mechanics can be formulated. This is supposed to be equivalent not to the 
bosonic c = 1 string, but to a worldsheet supersymmetric version, the c = 1 type OB 
string^, ||. We will not discuss this latter theory here. But many of the observations 
in the present review presumably can, and should, be generalised to the noncritical 
type OB (and OA) string background. 

In its simplest (un-orbifolded) form, the c = 1 string has a translationally invari- 
ant space/time direction X. If it is Euclidean, it can be compactified on a circle of 
radius R. In that case, all physical quantities (partition function and amplitudes) de- 
pend on R. The Euclidean theory can be interpreted as a finite-temperature theory 
with R labelling the inverse temperature. 

The value R = 1 (in units where a' = 1) is special because the c = 1 CFT is 
then self-dual under 

(i.D 

In this article we will focus on this background: the Euclidean c = 1 string with 
the X direction compactified at the self-dual radius. For short, we will refer to the 
theory as "c = 1,R = 1". 

There are many indications that string theory is topological in this background. 
The term "topological string theory" is usually taken to mean a string theory where 
the matter sector has a twisted M = 2 superconformal algebra with central charge 
zero. To make it a string theory, this matter then has to be coupled to topological 



worldsheet gravity [14], |15|| , which can be described by a set of bosonic and fermionic 
ghosts with total central charge zero. The fermionic charge of the twisted Af = 2 
algebra is the BRST charge defining physical states. 

A classic example of such a theory is a superstring on a Calabi-Yau background, 
whose superconformal worldsheet theory has Af = 2 supersymmetry with central 
charge c = 9. On making the standard topological twist of the superconformal 
algebra: 

T(z) - T(z) + \dJ{z) (1.2) 
where J(z) is the U(l) current, the algebra acquires central charge and we have a 



suitable matter system for a topological string fllS 



It is known that ordinary string theory, including the ghost sector, can be consid- 
ered as topological matter |L7|, |18|, [19[|, which makes the distinction between topologi- 
cal and non-topological theories less clear-cut. Even critical bosonic and superstrings 
are known to be "topological" in this sense. The best distinction one can make is 



-2- 



that the latter theories are "already twisted" while the conventional topological the- 
ories are formulated as TV = 2 superconformal theories and then given a topological 
twist. 

Topological string theories are generally related to the topology of the moduli 
space of Riemann surfaces. For example, "pure" topological gravity describes inter- 
section theory on cohomology classes associated to moduli space, or to vector bundles 
on moduli space [TH|. 

Some of the indications that c = 1, R — 1 is topological arise from its relation to 
other theories. The partition function of c = 1, R = 1 is closely related to that of the 



Penner matrix model |20|, |21|j , a model constructed to count the Euler characteristic 
of the moduli space of punctured Riemann surfaces. Amplitudes in c = 1, R = 1 
are summarised in the form of Woo constraints |22|] and the partition function of the 



perturbed theory is a r-function of an integrable hierarchy 1 . This r-function in turn 
can be written as a matrix model, the matrix model [El| (this has been previously 
referred to as the "Kontsevich-Penner model" [22[ and a "Kontsevich-type model for 



rg). 



Another indication of the topological nature of c = 1, R = 1 is that it is dual 



to the topological 2d black hole|T3, 25, EG]. The latter theory indeed starts life as 



an M = 2 superconformal field theory in two dimensions, where the second super- 
symmetry is a consequence of the Kazama-Suzuki construction. This theory is also 
dual to topological Landau- Ginsburg theory [^27], |28| , |29|, which is a convenient for- 



mulation for explicitly computing amplitudes and comparing them to those of the 
original c = 1, R = 1 theory test of the duality. 

Finally, c = 1, R — 1 has been shown to be dual to topological strings on a 
conifold singularity pOR. Recently it has also been argued that c = 1,R = oo, 



which might be more "physical" , is dual to an infinite-order orbifold of the conifold. 



Such infinite-order orbifolds can be understood in the language of deconstruction [p2 



33fl . In particular, the above orbifold is believed to be dual to the (2,0) CFT in 6 
dimensions ||34|| (an alternative limit instead gives rise to the nonrelativistic type IIA 
string |55|). Thus it may be that the noncompact c = 1 string can also be usefully 
formulated as a topological theory. In this case one could even try to continue the 
X direction back to Minkowski signature and see what happens in the dual theory. 
Progress on these issues is, however, quite limited to date. 

The various developments described above are important because the "topolog- 
ical" context is usually simpler than the "dynamical" one. Moreover, the former 
is often embedded in the latter, an example being Gopakumar-Vafa duality |36 



in 



which an open-closed topological string duality can be lifted |37j to an open-closed 
string duality in the full superstring theory. The computations in the topological 
theory are related to specific types of amplitudes in the dynamical theory. 

1 A different approach to constraints and associated matrix models is described for example 
in Ref.il. 
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In this article we will focus on the Penner and matrix models. These have 
been dubbed "topological matrix models" because they originate in the description of 
the moduli space of Riemann surfaces and also because they are dual to c = 1, R — 1. 
These are models of constant matrices and are quite different, in spirit as well as in 
details, from the familiar c = 1 matrix quantum mechanics. In particular, the nature 
and role of the large- iV limit is rather different. 

The plan of this article is as follows. In Section ^ we start by briefly reviewing 
the Liouville and the matrix quantum mechanics approaches to c = 1 string theory. 
After that we collect some explicit results from matrix quantum mechanics that will 
be useful in the subsequent discussions. 

In Section [3] we will explain how the Penner model is related to the triangulation 
of the moduli space of Riemann surfaces, rather than the triangulation of random 
surfaces. This requires some discussion of meromorphic quadratic differentials on 
Riemann surfaces. Section ^ deals with the model that is built from correlators 
of c = 1,R = 1. The partition function is a r-function of the Toda hierarchy. 
Relationships between this model and the Penner model as well as the Kontsevich 
model are discussed. Essentially the model is a perturbation of the Penner 
model by tachyon couplings, but we will see that the relation between them involves 
a somewhat subtle change of variables that may have a physical significance. 

To date, both these models have been constructed a posteriori by computing 
appropriate amplitudes of the c = 1, R — 1 string from matrix quantum mechanics, 
and then summing them up into an appropriate generating function. This is a 
somewhat ad hoc state of affairs and does not provide a fundamental explanation of 
how they arise, nor does it indicate how these models can be generalised to other 
backgrounds. 

In Section |5|, we will present some observations which might help to remedy this 
situation. It will be shown that the model can be rewritten as a "Liouville matrix 
model" - a model of constant matrices with an exponential plus linear potential. This 
way of rewriting the model is possible only because of its specific form, including the 
dependence on coupling constants. 

The presence of a matrix with an exponential plus linear potential is strongly 
suggestive of D-instantons moving in the background of a linear dilaton as well a 
cosmological (exponential) term arising from the closed string tachyon. We will 
try to develop the analogy, but will not give a conclusive argument that the Woo 
model is really the world-point action on iV D-instantons in the c = 1, R — 1 string 
theory. If this connection can really be demonstrated, it would indicate a new type 
of holographic relationship, analogous to the IKKT matrix model|| of critical string 
theory. It might also point to some new relationships between the perturbation series 
of string theory and the topology of moduli space. 

An important recent development is the construction of a new matrix model 
called the "normal matrix model" (NMM) . Like the model, this one too is 
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built directly out of the solution for the partition function of c = 1, R — 1 (and 
also other integer R). However, it does not seem to be known at present what is 
the precise relation of NMM to the and Penner models. We will briefly discuss 
the NMM at the end in order to exhibit some similarities to the other models that 
are the main subject of this article. Interestingly, in Ref.|38|] a duality has been 
proposed between matrix quantum mechanics and NMM, where the two models are 
associated respectively to the non-compact and compact cycles of a certain complex 
curve. This suggests an elegant answer to the question of how c = 1 matrix quantum 
mechanics and topological matrix models are related. The discussion of the normal 
matrix model presented here will unfortunately be very brief. 



2. c = 1 at R = 1 

In this section we will summarise some results from the matrix quantum mechanics 
approach to c = 1 string theory. An excellent review is Ref.|39|], where derivations 



of the formulae to be presented below, as well as extensive references, can be found. 
First we will briefly discuss the definition of c = 1 string theory, both from the 
continuum perspective and through matrices. 

2.1 Continuum formulation of c = 1 

The continuum description of noncritical c = 1 strings starts from a matter conformal 
field theory of a free scalar field X, with energy-momentum tensor: 

Tf z = -dXdX (2.1) 

This system has unit central charge, and its basic conformal fields are the vertex 
operators 

V = : e ikX : (2.2) 

as well as polynomials in d n X for various n, and products of these with the vertex 
operators. The coordinate X is taken to be the time direction, though one also 
frequently studies a Euclideanised theory where X is taken to be spacelike. 

A string theory with the critical central charge arises by coupling the above CFT 
to a Liouville field with energy-momentum tensor: 

= -d<j)d(f) + 2d 2 <j) (2.3) 

This Liouville field has central charge c = 25, so that the string is in this sense 
critical. However, it is supposed that this field arises as the scale factor of the 
worldsheet metric, which adjusts itself self-consistently to carry this central charge. 
The Liouville field is spacelike. 
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Observables of this theory are obtained by computing the BRST cohomology, 
just as is done for critical string theory. An important class of observables are the 
"tachyons" , which in Euclidean signature are given by: 

Tk = e (2-\k\)4> e ikx ( 2 A j 

The name "tachyon" is a misnomer, as these are actually massless states in two 
dimensions. 

There are also other physical modes in the BRST cohomology, the so-called 
"discrete states" POL ETl |3, which arise at special values of the momentum. These 



can be thought of as the two dimensional analogues of the graviton and other tensor 
modes of a closed string theory. In two dimensions, tensor fields have no propagating 
field-theoretic degrees of freedom, but in this theory they do have residual discrete 
modes that survive in the BRST cohomology. The tachyons described above can 
scatter into these discrete states, a fact which will become important when we discuss 
amplitudes. 

2.2 Matrix formulation of c = 1 

A very powerful description of c = 1 strings can be obtained starting from a random 
matrix integral: 



(2.5) 



where the matrix measure is given by: 

[dM] = Yl dM u Yl dMij dM*j (2.6) 

i i<j 

A perturbative expansion of this integral leads to 't Hooft-type Feynman diagrams 
with cubic vertices. Each such diagram specifies a unique Riemann surface topology 
on which it can be drawn as a random lattice with three lines meeting at every point. 
The faces are arbitrary n-gons. The dual lattice has n lines meeting at a point but 
the faces are triangles. The result is a triangulated Riemann surface. 

Thus the matrix quantum mechanics can be taken to represent the discretisation 
of a Riemann surface. Summing over all triangulations of the Riemann surface 
amounts to summing over all inequivalent conformal classes, so the result should be 
equivalent to a string theory. Indeed, this is the c = 1 string. The time direction 
that we called X in the continuum case is associated to the time coordinate t on 
which the matrix coordinate depends, while the Liouville direction in the continuum 



language is related to the infinite number of eigenvalues of the matrix at large N [33 

To solve this string theory, the matrix is first diagonalised. Then the eigenvalue- 
dependent part of the measure becomes: 

[dM} = Y[d\ l A 2 {\) (2.7) 
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Figure 1: A Riemann surface, with solid lines representing the matrix model diagram and 
dotted lines representing the dual triangulation. 

where A(A) = rL<j(-^ ~~ ^j) * s ^ ne Vandermonde determinant. In terms of eigenval- 
ues, the kinetic term of the Hamiltonian is: 

A^?i A(Ai) (2 - 8) 

acting on wave functions x(Aj). By redefining the wave function as follows: 

X (Xi) - *(Ai) = A(A l )x(A J ) (2.9) 

we have a simpler Hamiltonian but the new wave function \& is fermionic, since 
interchange of any two eigenvalues gives a minus sign. So the matrix eigenvalues 
behave as N fermions moving in the given potential. 

The ground state of the theory is found by filling up the first iV levels (after 
regulating the potential) to the Fermi level —/j,f- Next we take a double-scaling 
limit 

fi F ^0, f3-+oo (2.10) 

with /i — P/ip fixed. In this limit, for purposes of perturbation theory the eigenvalues 
behave like fermions in an inverted harmonic oscillator potential. 

The physical quantities (partition function, correlation functions) of this model 
will be given perturbatively as an expansion in powers of 




(2.11) 



Nonperturbative effects will typically arise as terms like 

e -l/9s„ e -tJ> (2.12) 

in the amplitudes. 

Computing the density of eigenvalues, one is able to evaluate the partition func- 
tion and free energy in a genus expansion. Correlators of observables can also be 
computed, as we will indicate in a subsequent section. 

This much of the discussion is valid for any compactification radius R, including 
the limit R — > oo. But at it! = 1 we will see that the formulae acquire special 
properties. 
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2.3 Free Energy 

Because computations in matrix quantum mechanics are highly technical, we will 
only quote some relevant results and refer the reader to Ref . [ 3J| and references therein 
for the derivations. 

The free energy = log Z{fi) of the c = 1 string was first obtained for 

arbitrary R by Gross and KlebanovfHJ: 



^-hf^.^- (2 . 13) 



dfi 2 J t sinht/2 smht/2R 

This is to be understood as an asymptotic expansion in powers of l/fi 2 = g 2 . The 
dependence on the compactification radius R arises from the last factor in the inte- 
grand. In the limit R — > oo, this factor goes to 1. 

The above formula can be integrated twice in fi to find the free energy. The 
integration constants are non- universal terms that will be unimportant. 

Performing the expansion in l/fi 2 , we find: 



fg(R) 
where 



fi~ 29 (2.14) 



k=0 



f g (R) = (2g - 1)! \2 2k - 2||2 2 ^) - 2^ ^^ R°~ 2k (2.15) 



and Bik are the Bernoulli numbers. 

Something special happens at R = 1. One can use a bilinear identity on Bernoulli 
numbers, due to Gosper 2 : 

" (1 - 2i-*) (1 - T- n + 1 )B n _ i B i _ (1 - n)B n 
^— ' (n — n\ 

Using this identity, we get: 

^ = _i og „ + f;^i| B2s |^ (2 . 17) 

r=i 2q 

g=l a 



dp 

Integrating, the free energy is found to be 



r{fi) R=1 = y 2gl a 2 - 29 (2.i8) 

V) j^2g{2g-2)^ 1 ] 



2 For just about everything you need to know about Bernoulli numbers, including this identity, 
see Ref.Ii. 
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For genus g — 0, 1 the coefficients are formally divergent (because of the dropped 
integration constants). 

It is remarkable that the above expression arises from doing matrix quantum 
mechanics in the double scaling limit. As we will see, it is closely related to the 
(virtual) Euler characteristic of the moduli space of Riemann surfaces. This is the 
first of many special properties of c = 1, R = 1 that we will encounter. 

2.4 Tachyon correlators 

Let us first consider the zero-momentum tachyon. The corresponding operator T = 
e 2 ^ is called the cosmological operator, because it is conjugate to the cosmological 
constant /i: 

(ToT ---T ) g ~ (2-19) 

Next consider T k for general nonzero k. At R = oo, A; is a continuous real variable, 
while at R = 1, it is quantised to be a (positive or negative) integer. Because A is a 
free field, the total k is conserved: 

n 

(T kl T k2 ■ ■ ■ T kn ) g = unless J* = (2.20) 

i=l 

In matrix quantum mechanics, one must identify the operators that correspond to 
the massless tachyon modes. This procedure is less straightforward than computing 
the BRST cohomology in the continuum case. It turns out that the correct local 
operators are defined as moments of loop operators^, |47|| , up to normalisation - 
about which we will say more below. Define: 

0(l,k) = J dxe ikx tr e~ m{x) ~ l k P k + ■ ■ ■ (2.21) 

where we keep the leading term for small I. 

Correlators of the P k the follow by computing correlation functions of the eigen- 



value density. For general R, these computations were performed in Refs. pq, ^2 
For example, the three point function is given by: 



(P kl P k2 P k3 ) = R5(h + k 2 + k 3 ) ^(IW 2 |+l*3|-2) r(1 _ N)r(1 _ N)r(1 _ \ h 



x 



^L(N - i)(N - 2) U(kl + k\ - n - 1) - \\ //-- i 



(2.22) 



It is a remarkable achievement of matrix quantum mechanics that this result is 
actually known to all orders in l//i 2 . Despite appearances, it is totally symmetric 
under interchange of ki,k 2 ,k 3 , and also invariant under fcj — ► — ki as required by 
translational symmetry in the A coordinate. 
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By contrast, amplitude computations in the continuum Liouville theory have 
only been done for genus g = ||49|| . The three-point function in this formulation was 
found to be: 



(T kl T k2 T k3 ) g=0 = R5(k l + k 2 + k 3 ) \x 



:(| fcl | + fel + fel-2) r(l-|feil)r(l-|A: 2 |)r(l-|A;3|) 

r(N) r(N) r(N) 

(2.23) 

It follows that the continuum tachyons T k and the matrix model puncture operators 
P k are related by: 

n = r^-p. (2.24) 

At zero momentum this represents an infinite change of scale between the two de- 
scriptions 3 . For other values of k it is just a finite change. 

We see that whenever A; is a nonzero integer, the amplitudes for both T k and P k 
develop poles. Integer momenta are precisely the ones for which additional discrete 
states (to which we briefly referred in the Introduction) exist in the theory. The 
divergences in tachyon amplitudes are believed to be due to the production of these 
states in intermediate channels. 

At R = 1, since k is always an integer, every correlator of T k or P k (except for 
Pq) is divergent. Hence it is necessary to define "amputated tachyon operators" T k : 

r(l + |*[) k P 9f(\ 

Tk - m-\k\) Tk - m-\k\) Pk (2 - 26) 

The correlators of these operators have no poles 4 . Indeed, the three-point function 
of amputated tachyons is: 

(T kl T k2 T k3 ) = 5(h + k 2 + k 3 ) ^\+\Ws\-2) \ klk2h \ x 

1 



24 (\k 3 \ - 1)(N - 2){k\ + k\ - |A; 3 | - 2) /i- 2 + 0(/i" 4 ) 



(2.27) 



The absence of poles in the amplitudes for amputated tachyons suggests that we 
are really working with a somewhat different theory, in which the other discrete states 
are no longer present. The discrete tachyons appear to be the only surviving states. 
In particular, this means that the radius perturbation, which should be turned on 
to continuously deform R away from its self-dual value, is absent. In this sense, the 



3 This is responsible for the fact that the cosmological operator is not really e 2 ^ but <\>e?^ . Indeed, 
we see that 

P Q = lim </>e 20 + infinite term (2.25) 

k^Q k 

4 In the Minkowski theory these leg pole factors are phases. Their physical role has been discussed 
by Polchinskipfl 
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theory is "stuck" at R = l 5 The above properties characterise the topological c = 1 
background 6 . 

It is intriguing that the coefficient of the fi~ 2 term is an integer for every k±, &2, ^3 
satisfying k\ + k 2 + k 3 = 0. Explicit computation of higher-order terms in this and 
other correlators suggests that this is always the case, though there does not seem to 
be an explicit statement or proof in the literature. If every term is an integer, this 
would certainly be special to R = 1, as this property can be seen not to hold for any 
other R. 

An exact generating function for all discrete-tachyon correlators, in every genus, 



is known |22|. This will be described in a subsequent section. 



3. Riemann Surfaces and the Penner Matrix Model 

As has been indicated, the above results were derived starting from matrix quantum 
mechanics, using the free fermion description. The first indication that these are 
related to the topology of Riemann surfaces arose in the study of the Penner matrix 
model|2lJ, via the work of Distler and Vafa[|21||. To explain this model, we must first 



discuss some mathematical issues relating to the moduli space of Riemann surfaces 
of genus g. 

3.1 Moduli space of Riemann surfaces and its topology 

The moduli space of a compact Riemann surface of genus g without punctures, -M 9 , 
is a manifold of complex dimension 3g — 3 (for g > 1). It arises as the quotient of 
the Teichmuller space by the action of the mapping class group. Since this action 
has fixed points, the moduli space A4 g has orbifold-like singularities, 

What is the simplest topological invariant of For a D-dimensional smooth 
manifold, we can define the the Euler characteristic x by making a simplicial decom- 
position S of the manifold and then evaluating: 

X = £(-1)* 

where di is the dimension of the ith simplex, and the sum is over all the simplices in 
the complex S. 

In the presence of orbifold singularities, the natural quantity to define is the 
virtual Euler characteristic xv- Here each term in the sum is divided by the order 
of a discrete group Tj that fixes the ith simplex. Thus: 

xv = S (3 - 2) 



5 One can take orbifolds to get other integer values of R. 

6 Note that "winding tachyons" , which are winding modes around the X direction, should still 
exist, though they have not yet been satisfactorily understood in topological matrix models. 
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It was found by Harer and Zagier[[5l|1 that for the moduli space of unpunctured 



Riemann surfaces, the virtual Euler characteristic is: 



* = (3 - 3) 

where B 2g are the Bernoulli numbers. 

One can also consider the moduli space M. 9lS of Riemann surfaces of genus g 
with s punctures. This has complex dimension 3g — 3 + s (for g>l,org = l,s>l, 
or g = 0, s > 3). Its virtual Euler characteristic is given by: 

(-l)'(2g-3 + 8 )!(2g-l) _ (2g-3 + s\ 

X 9 ,s- S^-(-l) ^ s jx, (3.4) 

The above results for \g an d Xg,s are related in a suggestive way. In string theory, 
one expects that s punctures can be created by differentiating the vacuum amplitude 
s times with respect to the cosmological constant. And indeed, as noticed by Distler 



and Vafa||21||, one has the identity: 

\^s(x 9 ^ 2a )=X 9 ,s^- s (3.5) 

Thus, if we could invent a model whose genus g partition function is \gi we might 
expect it to bear some relation to string theory. Below we will see that this is indeed 
the case. 

3.2 Quadratic differentials and fatgraphs 

The above results were obtained by triangulating the moduli space of punctured 



Riemann surfaces. Such a triangulation was constructed by Harer in terms of 
quadratic differentials, using a theorem due to Strebel|53[. It is instructive to sketch 
how this was done. 

On any Riemann surface with a finite number of marked points, one can define 
a meromorphic quadratic differential 

V = Vzz(z) dz 2 (3.6) 

which has poles at the locations of the marked points. This can equivalently be 
considered as a holomorphic differential on the punctured Riemann surface, which 
has these marked points removed. 
Under a change of coordinates: 

z -> z'(z) (3.7) 
a quadratic differential transforms as: 



f dz\ 



2 



V'z'z'(z') = J Vzz{z) (3. 
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This differential can be used to invariantly define the length of a curve 7 on the 
Riemann surface. The length is simply defined to be: 



7l 



n 



V\vW\\dz\ (3.9) 



7 



Indeed, defining a new coordinate via 

dw = \fq{z)dz (3.10) 

we see that this length is the ordinary length of the curve in the Euclidean sense, in 
the w coordinate. 

Now a quadratic differential is not in general unique, but it was shown by 
Strebel|)3|] that with some extra properties, such a differential is unique upto multi- 
plication by a positive real number. For this, let us consider a geodesic curve under 
the metric defined above. At any point, such a curve will be called horizontal if 77 
is real and positive along it, and vertical if 77 is real and negative. The horizontal 
curves define flows along the Riemann surface. 

The flow pattern is regular except at zeroes and poles of 77. Here the flows exhibit 
interesting properties. At an rath-order zero of the quadratic differential, precisely 
n + 2 horizontal curves meet at a point. To see this, let us write down what the 
differential looks like near this zero and along the radial direction: 

77 ~ z n (dz) 2 ~ e i{n+2)e dr 2 (3.11) 

From this we see that as we encircle the zero, there are precisely n + 2 values of the 
angle 9 at which this differential is positive, or horizontal. 

At a double pole, if the coefficient is real and negative, the flows form concentric 
circles around the point. This follows from the fact that near the pole, and along 
the angular direction, the differential looks like: 

d 7^ 

v ~ - c —~cd6 2 (3.12) 

z 2 

Thus, in the 9 direction, the differential is positive, or horizontal, at all points sur- 
rounding the double pole. 

Other behaviours are possible at poles of order n = 1 or n > 3, or if the coefficient 
of 77 at a double pole is complex (for a comprehensive description, see Figs. 5-13 of 
Ref. f53[). But we do not need to consider these other behaviours because we will 
restrict our quadratic differentials not to have such poles. 

In fact, we will specialise to a class of quadratic differentials which have a double 
pole at a point P, at which the coefficient c is required to be real and negative. 
Next, we require that all smooth horizontal trajectories (i.e., those that do not pass 
through zeroes of 77) form closed curves. Quadratic differentials satisfying all these 
conditions exist, and are called horocyclic. 
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An illustration of the flow pattern associated to a horocyclic quadratic differential 
is given in Fig.2. Note that the vertex has five lines meeting at a point, indicating a 
third-order zero (the differential also has additional zeroes not shown in the diagram). 



Double pole with negative coefficient 




Third order zero 



Figure 2: Riemann surface with the flow pattern of a horocyclic quadratic differential. 



Strebel's theorem states that on every Riemann surface of genus g with 1 punc- 
ture, for fixed values of its moduli, there exists a unique horocyclic quadratic differ- 
ential with a double pole at the puncture. (The uniqueness is upto multiplication by 
a real positive number). 

Thus, by studying how these quadratic differentials vary as we vary the moduli, 
we get information about the moduli space M. g ,i of a once-punctured Riemann sur- 
face. In view of some mathematical results relating the moduli spaces of Riemann 
surfaces of a fixed genus g but different numbers of punctures sfl52"l, this construction 
actually enables us to study the moduli space Ai g , s of s-punctured Riemann surfaces 
as well. 

The main motivation for having gone through this mathematical description here 
is that we will now see the emergence of "fatgraphs", otherwise known as 't Hooft 
diagrams or large- N diagrams. Most of the flows are closed and smooth, but there 
are singular ones that branch into n + 2-point vertices at rath order zeroes of r\. We 
can think of these singular flows as defining a Feynman diagram, whose vertices are 
the branch points. Each double pole of rj is a point around which the flows form 
a loop, hence the number of loops of the diagram is the number of double poles, 
which is the number of punctures of the original Riemann surface. Finally, because 
the flows that do not pass through a zero are smooth, each singular flow can be 
"thickened" into a smooth ribbon in a unique way, and we arrive at a fatgraph. 

The fatgraphs with a single face triangulate the moduli space M g ,i in the follow- 
ing way. Consider the lengths of each edge of a fatgraph, as computed in the metric 
defined in Eq. (|3.9| ) above. We have an overall freedom of scaling the whole Riemann 
surface, which clearly does not change the moduli. So to vary the moduli, we must 
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change the lengths of the different edges keeping the total length fixed. By Strebel's 
correspondence between quadratic differentials and moduli, this sweeps out a region 
of the moduli space of the Riemann surface. The dimensionality of this region will 
be E — 1 where E is the number of edges of the graph. 

As Harer argues, this region swept out will not be the whole moduli space, but 
a simplex of it. In a simplicial decomposition, at the boundary of a simplex we find 
a lower-dimensional simplex. In terms of fat graphs, a boundary occurs whenever 
a length goes to zero and two vertices meet. As an example, when two three-point 
vertices meet, we obtain a four-point vertex. The new graph thus obtained sweeps 
out a lower-dimensional simplex in the moduli space as its remaining lengths are 
varied, and this process continues. 

Now the virtual Euler characteristic can be defined directly in terms of fatgraphs. 
We consider the set of all fatgraphs of a given genus g and a single puncture. Let 
us call the set S and label each distinct graph by an integer i e S. Let Tj be the 
automorphism group of a fat graph. Then, defining dj = E — 1, we write 

which is analogous to Eq. (|3.2| ) above, except that now the sum is over fatgraphs 
rather than over simplices. In particular, the automorphism group of the fatgraph is 
the same as the group that fixes the corresponding simplex. Very low-dimensional 
simplices do not contribute to xv since at some point the corresponding fatgraph has 
an automorphism group of infinite order, in which case the denominator of Eq. ( |3.13| ) 
becomes infinite and the corresponding term vanishes. 

Let us check how the correspondence between fatgraphs and quadratic differen- 
tials works out in practice. The fatgraphs we have been considering have V vertices, 
E edges and 1 face. These integers satisfy: 

V-E + l = 2-2g (3.14) 

where g is the genus of the Riemann surface on which the graph is drawn. 
We also have the relations 

v = Y,v k , E = ^j2 kv * ( 3 - 15 ) 

k k 

where is the number of fc-point vertices. From these relations, we get: 

J2(k - 2)v k = Ag - 2 (3.16) 

k 

All integer solutions of this equation, i.e. all choices of the set {vk\ for fixed g, are 
valid graphs that correspond to simplices in the triangulation of M. 9j x- 
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Let us recast the above equation as 



J2(k-2)v k -2 = 4g-4 (3.17) 

k 

Then, noticing that A; — 2 is the order of the zero represented by a Appoint vertex, 
we see that the first term on the left is the total number of zeroes (weighted with 
their multiplicity) of the quadratic differential corresponding to the given fatgraph. 
Moreover, the differential has precisely one double pole, so the second term is minus 
the (weighted) number of poles. Thus this result is in accord with a theorem which 
states that for meromorphic quadratic differentials on a Riemann surface of genus g, 

#(zeroes) - #(poles) = Ag - A (3.18) 

A particular solution that is always available is v$ — V, — 0, k > 4. Clearly 
this gives the maximum possible number of vertices and therefore also edges. In this 
case, 

V = 4g-2 (3.19) 



and we see from Eq. (|3.15 ) that the number of edges is 



E = ^V = 6g-3 (3.20) 

Thus the dimension of the space spanned by varying the lengths of the graph keeping 
s lengths fixed, is: 

E- 1=60-4 (3.21) 

which is the real dimension of A4 g ,i- Thus the graphs with only cubic vertices span 
a top-dimensional simplex in moduli space. All other graphs arise by collapse of one 
or more lines, merging two or more 3-point vertices to create higher n-point vertices. 
These correspond to simplices of lower dimension in the moduli space. 

As an example, in genus 3 with one puncture, the moduli space has real dimension 
14. Hence the graph that describes the "bulk" of moduli space should have 15 edges 
and 10 3-point vertices. Also, because there is a single puncture, the graph should 
have just one face. In summary, we have 

V = 10, E = 15, s = 1 V-E + s = -A = 2-2g (3.22) 

as expected. 

3.3 The Penner model 

In 1986, Penner [pOfl constructed a matrix model that provides a generating functional 
for Xg,s- The Penner model is defined in terms oi N x N random matrices whose 
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"fatgraphs" are precisely the ones described in the previous subsection. The free 
energy T = log Z of this model then has the expansion: 

where t is a parameter of the model. The term s = is not present in the sum. 

Let us first present the model and then show why it is correct. The model is 
given by an integral over Hermitian random matrices: 

Z = M P J [dQ] e~ Nt tr ^"=2 i ® k (3.24) 

= Mp J [dQ] e Nt tr ( Ml ~Q)+Q) (3.25) 

where Mp is a normalisation factor given by: 

Up 1 = J [dQ] e~ Nt tr \ Q2 (3.26) 

and the matrix measure [dQ] is given by: 

[dQ] = JJ dQ ti J] dQa dQ% (3.27) 

i i<j 

This action has all powers k > 2 of the random matrix appearing in it. The model 
is to be considered as a perturbation series around Q ~ 0. 

To show that this model is correct, we must show that its fatgraphs are in one 
to one correspondence with those arising from quadratic differentials as discussed in 
the previous section. Thus we must show that the partition function involves a sum 
over all (connected and disconnected) fatgraphs for a fixed genus g and number of 
faces s, multiplied by the weighting factor 

N 2~2 H 2-2 9 -s = { _ Nt) V m -E {Nr {328) 

By taking the logarithm of the partition function we end up with the sum over 
connected graphs in the familiar way. 

In the above expression Tj, the automorphism group, is the collection of maps 
of a given fatgraph to itself such that (i) the set of vertices is mapped onto itself, 
(ii) the set of edges is mapped to itself, and (iii) the cyclic ordering of each vertex is 
preserved. A key result due to Penner [p0|1 is that the order of this group is equal to: 

= U l* c ( 3 - 29 ) 

vertices 

where C is the combinatoric factor labelling how many distinct contractions lead to 
the same graph. Now this is exactly the factor that arises if we obtain our diagrams 
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as the expansion of the matrix integral Eq.( [3.24 ), with the V\ coming from the order 



of expansion of the exponent, the \ factors being built into the action, and the 
combinatoric factor being the usual one. 

The remaining factors are straightforward. Each vertex has a factor —Nt, each 
edge is a propagator given by (Nt) -1 , and each face gives the usual index sum N. 
This completes our informal derivation of the Penner model 7 . 

3.4 Double-scaled Penner model 

In principle, the Penner model has nothing to say about the moduli spaces M. g fl of 
unpunctured Riemann surfaces. As we have noted, the expansion of the partition 
function does not contain terms with s = 0. 

However, it was noticed by Distler and Vafa [plJ| that in a suitable double-scaling 



limit, the model seems to describe unpunctured Riemann surfaces. Their observation 
goes as follows. Start with the equations: 

T = T a = £ %m N2 ~ 29 t2 ~ 29 ~ S (3-30) 

9 9,s 

(-l)'(2g-3 + a)!(2g-l) p [2g-3 + s\ 

X 9 ,s = B 2g = (-1) ^ s j Xg (3.31) 

where we have inserted the known value of Xg,s- 

Let us work with genus g > 1. It is straightforward to explicitly perform the 
sum over s, and one gets: 

2g(2g - 2) 

Now Distler and Vafa took the limit iV — > oo and t — > t c = — 1, keeping fixed the 
product N(l +t) — —v. This leads immediately to the result: 



^^2^2T^ 29 (( 1 + ^- 1 ) ^ 



T a = B2g - s v 2 - 29 = Y X y- 29 (3.33) 
9 2g(2g-2) ^ Xg K } 

and we have recovered the virtual Euler characteristic of Ai g .o ] - 

Thus the Penner model, originally designed to study the moduli space of punc- 
tured Riemann surfaces, describes unpunctured ones too - but only in the special 
double-scaling limit above. We can remove some of the mystery of the double-scaling 
limit by defining, for any finite N and t, the parameter 

v = -N-Nt (3.34) 

in terms of which the matrix model can be written: 

Z = M P J \dQ] e^ +N) tr ( ~Q)+Q) (3.35) 



7 Another derivation, as well as proposed generalisations, can be found in Ref.|54 
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Now the double-scaling limit is simply the ordinary large- N limit, N — > oo, with v 
held fixed. 

From the discussions of the previous sections, it should be evident that the 
double-scaled Penner model is closely related to c = 1, R — 1. Its free energy: 

^E^r 2 - 29 < 3 - 36 > 

is almost identical to the free energy of the c = 1, R = 1 string theory: 

^^dtV" < 3 - 37 > 

However, there is an important issue of signs. It is well-known that the Bernoulli 
numbers alternate in sign: 

\B 2g \ = (-ly-'Big (3.38) 
Therefore if we define v = —ifi, we can write: 

°° I R I 

^p(-iAi) = g 2 ^ 2) /i 2 " 29 = ^ c= i(/i) (3.39) 

Thus the double-scaled Penner model (at imaginary cosmological constant) has the 
same free energy as the c = 1, R — 1 string. This intriguing correspondence calls for 
an explanation. 

In the last section we will try to argue that the double-scaled Penner model could 
be related to the Euclidean c = 1, R = 1 theory, and could perhaps be derived from 
it, as the matrix integral for iV D-instantons of the latter theory. If true, this might 
both explain the above coincidence and also provide a role for D-instantons of the 
noncritical c = 1 string. 



4. The Matrix Model 



Let us now return to the solution of the c = 1,R — 1 matrix quantum mechanics. 
We will now focus on its tachyon amplitudes. 

4.1 c = 1 amplitudes and 

As was already mentioned, for this theory a generating functional for all tachyon 
correlators to all genus has been obtained^]. This functional F{t, t) depends on an 
infinite set of couplings tf-,ik such that: 



(4.1) 

t=t=o 
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where on the LHS we have connected amplitudes. We see that the couplings tk are 
sources for tachyons of positive momentum k, while the couplings tk are sources for 
tachyons of negative momentum —k. These tachyons are "amputated" in the sense 
of Eq.(CT). 

Instead of constructing the generating functional directly in terms of tk,ik, it 
turns out necessary to encode the parameters tk into a constant N x N matrix A. 
The tk are defined in terms of A by the relation: 



U 



vk 



trA~ 



(4.2) 



which is sometimes called the Kontsevich-Miwa transform. As iV — > oo, this matrix 
can encode infinitely many independent parameters tk- 

Now using matrix quantum mechanics at R — 1, it was shown [p2| that Z(t,t) = 
e F(t,t) j s a r-function of the Toda hierarchy, satisfying an infinite set of constraints 
that form a algebra. These constraints were subsequently rewritten in the fol- 



lowing form[p4H, which is the version that we will need for the subsequent discussion: 



1 dZ 

-u) dt n 



-v) 



■(det A)" tr 



d_ 

dA 



(det A)-"Z(t,i) 



(4.3) 



Here v = —ifi, and \i is the cosmological constant. The correlators determined by 
the above expression are actually invariant under hi — > —hi, though this is far from 
manifest, since t n and t n do not appear symmetrically. 



The Woo constraints can be integrated to give a matrix model |24] as follows . 
Let us start by assuming that Z(t,t) is an integral over constant matrices M: 



Z(t, t) = (det A) v J [dM] e tr ^( M ' A ' ^ 



(4.4) 



where [dM] = Yli dMu Yli<j dMijdM*j. The potential V is determined by imposing 
the above differential equation: 







_(-!/) 8t n 

This determines: 



tr 



d_ 

dA 



[dM]e trV(M,A,i) =Q 



oo 

V(M, t, i) = -v (MA + ^ Mfc ) + f( M ) 

fe=i 



(4.5) 



(4.6) 



where f(M) is a function independent of A, i. 

The function f(M) can be determined using a boundary condition, arising from 
conservation of the tachyon momentum. Momentum conservation tells us that if we 

8 The same constraints were integrated in Ref.[^2j but owing to some technical errors, as explained 
in Ref.p4], the resulting matrix model obtained in Ref.p2] was not correct. 
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set all the equal to zero, then Z(A, 0) has to be independent of A. From Eqns. 
and (14.61), we get: 



Z(t, 0) = (det A) u J [dM] e~ u tr MA + tr (4.7) 

Upon changing variables M — > MA' 1 , we see that: 

[dM] -»■ (det A)~ N [dM] (4.8) 

It follows that: 

Z(t,0) = (det AY~ N J[dM]e~ vtrM + tr /(MA -1 ) (4Q) 

\dM] e~ v trM + tr /(MA' 1 ) + (v - N) tr log A (41Q) 



This uniquely determines the function f(M) to be: 

f(M) = (v - N) log M (4.11) 

Putting everything together, we see that the generating function of all tachyon 
amplitudes in the c = 1, R — 1 string theory is: 

Z(t,t) = (det A)" J[dM]e ir ^ MA+ ^~ N ) l ^ M -^Z^M k ) {A12) 
= J [dM] e tr (~ uM + (v — N)logM — v Er=i 4(MA-!) fc ) (4 13) 

This is the matrix model [El]. The second form was obtained by redefining 
M -> MA' 1 . We also recall that 

t fc = i-trA- fc (4.14) 
v k 

Note that for the matrix integral to be well-defined, the matrix M (more pre- 
cisely, its eigenvalues) must be positive semi-definite. 

The Woo model has a number of interesting properties. One of these is the 
puncture equation, which says that the shift 

A^A-el, tx^ti + e (4.15) 

changes the free energy T by a simple additive factor: 

oo 

T(A - el, i k + 4,i e) = T(A, t) - v 2 ^ e k t k (4.16) 

fc=i 

This is easiest to see when the model is written in the form Eq. (|4.12| ). We see that 
there are two linear terms in M, one multiplying A and the other multiplying t\. 
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Their variations under the above transformation cancel out, leaving only the effect 
of varying the determinant factor in front of the integral. Taking logarithms gives 
the above behaviour of the free energy. 

Another property of the Woo model is that it is invariant under the simultaneous 
rescaling tk — > X~ k tk and tk — > X h ik, for any (in principle, complex) parameter A. This 
is just momentum conservation. To see this invariance, note first that by Eq. ( f4.2| , 
the rescaling — ► X~ k tk is the same as A — > A A Now the expression Eq.( 4.13|) is 
manifestly invariant under A — ► XA, ik — > X k ik- 

There have been other matrix models in the literature which have a form similar 
to that of the model 9 . Some examples are the Chekhov- Makeenko model [EC 



the NBI matrix model p7\, pq] and the normal matrix model|38|]. The first of these is 
similar in structure to our Woo matrix model but the coefficients are different. The 
second one is proposed as a modification of the IKKT matrix model by coupling it 
to a new matrix with a log plus linear interaction. The third example, the normal 



matrix model, will be discussed in Sec. |513. 



4.2 Relation to the Penner model 

The matrix integral of the Woo model is convergent for real positive v. Indeed, it is 
the matrix analogue of the T-function. But we are presently interested only in its 
perturbative expansion in inverse powers of v 2 . If we send v — > —u, we get the same 
perturbation series, even though the integral is formally no longer convergent and 
has to be defined by analytic continuation (like the T-function itself). 

It turns out that this transformation is necessary in order to recover the Penner 
model from the Woo model. Carrying out the change v — ► —u, the resulting matrix 
integral is: 

Z(t, 0) = J [dM] e tr ( uM ~(v + N ) lo § M ) (4.17) 
Now make the change of variables: 

M = a(l-Q) (4.18) 

where a is a parameter to be determined. 

Thus the W^ matrix model is transformed into the following matrix integral: 

Z = N J [dQ] e~( u + N ) tr Ml ~Q)-vaQ (4 ig) 

where M is a normalisation constant. 

Choosing a = 1 + N/ u, we get the standard form of the Penner model: 



Z=M J[dQ]e Nt tr(log(l-Q)+Q) 



(4.20) 



3 An early attempt to incorporate correlators in the Penner model is Ref. pi 
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where t — — 1 — vjN . 

Notice that in the "double-scaled" limit of the Penner model, we have iV — ► oo 
with v finite, so the parameter a goes to infinity. As a result, the Penner model 
perturbation series, around Q ~ 0, actually corresponds to the model in the 
region of large M. 

A word about normalisation is called for. The model determines correlators 
of tachyons of nonzero integer momentum (positive or negative). If insertions of 
Tq are also required, they can be obtained by differentiating in v. However, in the 
derivation leading up to Eq.( [4.1$ ), the overall ^-dependent normalisation was not 
fixed, as this does not affect correlators for which at least one tachyon has nonzero 
momentum (by momentum conservation, this means that at least two tachyons have 
nonzero momentum). 

This is why Eq.( |4.20| ) has a normalisation constant in front of it. We could of 
course have at the outset chosen the normalisation in Eq.( 4.13 ) to reproduce the 
correct normalisation factor of the Penner model given in Eq.( |3.26| ). But amusingly, 
this is not really necessary, for the following reason. We can extract correlators 
of zero-momentum tachyons from the model by taking a formal limit k — > 
in the corresponding combinatorial formulae. That this gives the right answer is 



demonstrated in Ref. 29 



4.3 Relation to the Kontsevich model 



We have not so far discussed the original Kontsevich matrix model|59|], but it is 
appropriate to make a brief mention of it here. Like Penner's model, this model 
too was formulated to solve a combinatoric problem, this time concerning stable 
cohomology classes on moduli space. Via Witten's results relating these cohomology 



classes to pure topological gravity and c < 1 strings [0], this matrix model actually 
describes c < 1 (rather than c = 1) noncritical strings. 

The Kontsevich model is described by random matrices with the following matrix 
integral: 

Z K {t) = J [dX] e~ tr l AX2 + % tr X3 /6 (4 . 21) 

The ti on which the partition function depends are an infinite set of parameters 
defined by: 

t k = -(2k - 1)!! trA" 2 ^ 1 (4.22) 

Differentiating the free energy Tk = log 2k with respect to the ti gives rise to an 
infinite set of correlation functions of operators that are analogous to the "amputated 
discrete tachyons" of previous sections. In turn, these correlators are the amplitudes 
of a class of c < 1 string theories based on (2, q) minimal models, the ones originally 
studied as double-scaled matrix modelsfBOl, BU, B!|. 
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It was noticed by Kontsevich that his model is related to the matrix Airy func- 
tion: 

i tr (X 3 /3 - XF N 



AiY) = J [dX] e V ' J (4.23) 

Indeed, the change of variables: 

X = 2- 1/3 X + Y 1/2 , Y = -2~ 2/3 A 2 (4.24) 

brings Eq.( |4.23| ) into the form of Eq.( |4.21|) , with, however, a significant prefactor 
depending on A. Thus the perturbation series of one can be expressed in terms of 
the perturbation series of the other. 

It was also noted in Ref.|SD| that the matrix Airy function Eq.( |4.23| ) is to be 



viewed as an asymptotic expansion at large Y, for which the saddle point at X ~ 
Y2 dominates and therefore the matrix X is also large. On the other hand, the 
corresponding perturbation series for the model of Eq.( [4.2l[) is around X ~ 0. This 



is strikingly similar to the relation between the model of Sec. |4], and the Penner 



model of Sec. |373|. As noted in Sec. |4^, the former model expanded about large 
values of the matrix M gets mapped to the latter model expanded about Q ~ 0. 
We can actually obtain the matrix Airy function directly as a special case of the 
model of Eq. (|4T3| ). Starting with Eq.( ^T3|) at finite N, set v = N and i k = 5 kj3 . 
The result is: 

Z(t) = (det Af IdM e~ N tr MA ~ N tr M " (4.25) 



After a further rescaling of both M and A, this is just proportional to the integral 
Eq.( f4.23| ). Thus the matrix Airy function arises from the model by "condensing" 
a specific negative- momentum tachyon, £3. 

Generalisations of the Kontsevich model to describe noncritical strings based on 
(p, q) minimal models do exist |63|]. They are related to generalisations of the matrix 



Airy integral where the cubic term is replaced by a higher power. Clearly these too 
can be obtained from the model by condensing, or turning on, ik for k > 3. 

For string theory amplitudes, we are interested not just in the matrix integral 
but also the A-dependent or A-dependent normalisations. So these observations do 
not show that c < 1 strings are obtained from the c = 1 background by "condensing" 
discrete tachyons, but they do seem to suggest this. Specifically, turning on i p or 
condensing T_ p is related to the family of (p — 1, q) minimal models coupled to 
gravity, for fixed p and all coprime q. More work is needed to precisely establish this 
relationship. 

5. Liouville Matrix Model and D-Instantons 

In this section we will rewrite the model in a form that involves a Liouville-like 
potential. This form will be suggestive of a D-instanton interpretation. 
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5.1 Unperturbed Woo and the Liouville Matrix Model 

Let us start with the model as written in Eg. (14.17 ), where all the have been set 



to zero. Notice that the coefficients in the action have a very specific dependence. 
The linear term is multiplied by a constant v that remains fixed in the large- N limit, 
while the log term has a coefficient v — N. Also, as we have seen, the matrix M has 
to be positive semi-definite. Both these features look slightly unnatural in a matrix 
model, and suggest a redefinition of variables. 
Define a new N x N matrix $ by: 

M = e $ (5.1) 

Under this transformation, the matrix measure transforms as: 

[dM] {detM)- N = (5.2) 

and the integral in Eq.( |4.17p becomes: 

Z = J [d$] e~ S = J [d*] e~ v tr ( e * - $ ) (5.3) 

where [d$] is an (analytically continued) matrix measure appropriate for unitary 
matrices: 2 

[d$] = J] dxf>i 11 (sinh htfi - <j>j)\ (5.4) 

i i<j 

Since M is Hermitian and positive semi-definite, it follows that $ is Hermitian with 
no further restrictions. Thus we no longer have to deal with a positivity constraint. 
More remarkably, we have also got rid of the explicit N in the exponent. The matrix 
model in the $ variable is simply a matrix integral with a potential that has a 
Liouville type exponential term, as well as a linear term. We will refer to it as the 
Liouville matrix model. 



5.2 Interpretation of the Model 

How do we interpret this model? The potential looks very much like that on the 
worldsheet of a fundamental string in the c = 1 background. If we think of <l> as the 
Liouville field, the two terms are reminiscent of the Liouville potential (cosmological 
operator) and the linear dilaton potential respectively. 

However, there is a different and more plausible interpretation. The coefficient of 
the entire matrix action is u, the cosmological constant. In noncritical string theories 
this corresponds to the inverse string coupling: v = l/g s . Thus the action of the 
Liouville matrix model has a factor l/g s in front of it. This is highly suggestive of a 
D-brane action. Since the matrices are constant, the D-brane in question should be 
a D-instanton. 
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In critical superstring theory, a single D-instanton has essentially no action cor- 
responding to it. But N D-instantons do have an action, consisting of all the com- 
mutator terms that arise in the nonabelian Yang-Mills action: 

S=-tr[A fl ,A u ] 2 + --- (5.5) 
9s 

Fermionic terms and higher-order commutators are also present. All these commu- 
tator terms arise because the scalars transverse to the instanton, and their fermionic 
partners, are promoted to matrices when there are N instantons. Systems of D- 
instantons for large N are associated to the IKKT|| matrix model of type IIB string 
theory. 

Suppose we had N D-instantons in c = 1 string theory. What action should 
we propose for them 10 ? One could try to address this by constructing a boundary 
state and computing open-string amplitudes. This approach would work only at 
weak coupling, that is far out near the region of large negative Liouville field <fi. In 
this region, the D-instanton will feel two forces, one due to the cosmological term, 
i.e. the Liouville potential, which is an excitation of the closed-string tachyon, and 
the other due to the linear dilaton. The former would drive the brane towards weak 
coupling, while the latter drives it to strong coupling. It is at least plausible that 
the two competing forces have an equilibrium position somewhere in the region of 
~ 0. However, it would be hard to gather reliable information in that region, since 
the theory is strongly coupled there. 

Hence we leave this question for future work. What is quite striking, though, 
is that without making any assumptions at all, we have obtained a Liouville matrix 
model which precisely consists of a competing exponential and linear potential. This 
might be considered as evidence for a D-instanton interpretation of the matrix 
model. 

This interpretation also adds something to the discussion in the previous section 
on the relationship of the Woo to the Penner model. Recall that the transformation 
from the to the Penner matrix, Eq.( [4.18D involves a large parameter a = 1+N/ v, 



such that at large N, small values of the Penner matrix Q ~ are mapped to large 
values of M. Now we see that in the Liouville variable this is the region of large <3>, 
or the strong coupling region. This could be related to the fact that in the DO-brane 



description of c = 1 strings[@] via the ZZ boundary state|L3]], the branes are thought 
of as being localised at large (p. 

Some ingredients are clearly missing from the story. We did not find a reason why 
the cosmological constant in this theory should be imaginary. Also, a D-instanton 
action should depend on two additional variables: the transverse scalar X describing 
the c = 1 direction, and the open-string tachyon T p e n- We will see that something 
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like X does appear when we consider tachyon perturbations, but we do not seem to 
find a role for T pen, unless it is somehow "mixed up" with $. 



5.3 Tachyon Perturbations 

We proposed that the matrix model on iV D-instantons is equivalent to the 
model without tachyon perturbations (i.e. at t,i — 0). What happens when we turn 
on tachyon perturbations? 

If the proposal is correct, the D-instanton action should match the full action 
as a function of t, t: 

Z(t, t) = J [dM] e tr (-^ M + ( u ~ N)\og M — v £~ > h{MA^f) (56) 

where 

t k = \trA~ k (5.7) 
v k 

Let us focus on the matrix A. If it is Hermitian then (for real v) the are real. 
But we could equally well choose A to be unitary, in which case the t k are complex, 
since we only need to differentiate jF(t, i) near t, t — 0. 

Let us therefore choose A to be unitary, and parametrise it as: 

A = e lX (5.8) 

where X is Hermitian. 

Now, in the Liouville-like variable the matrix model with tachyon per- 
turbations can be written: 



r c r -utr ( e *-<S> + Y°° t k (e*e- iX ) h ) 



(5.9) 



with 

t k = — tr(e~* kX ) (5.10) 
vk 

To understand this better, let us consider the simplest case of N — 1, i.e. lxl 
matrices. In this case, the perturbing term becomes: 

oo oo 

( e ^) fc ->5>e*<*-»> (5.11) 
fc=i fc=i 

The RHS looks like half of the mode expansion of a 2D Euclidean scalar field. This 
suggests that we identify x with the (compact) Euclidean time direction. Indeed, 
the periodicity is correct, since exp(ikx) is single-valued at R — 1 (i.e. x — > x + 27r) 
precisely for integer k. It is plausible that this is how a single D-instanton couples 
to closed-string tachyons of negative momentum. 
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For N D-instantons we need to go back to the noncommuting matrices, i.e. make 
the replacement: 

e K4>-ix) _^ tr ( e * e -^) fc (5. 12 ) 

This amounts to a specific prediction for the matrix ordering on N D-instantons of 
c = 1, R = 1 string theory. When the exponentials on the RHS are expanded, we 
find an infinite sequence of commutator terms between $ and X. 

The most puzzling aspect of this framework is that X is not a dynamical variable 
unlike $ (i.e. a random matrix to be integrated over), but rather a fixed background. 
It seems to be determined self-consistently, like a condensate, in terms of the positive 
momentum tachyon couplings, via the Kontsevich-Miwa relation Eq. ( |5.10| ). 

We can get an indication how this might come about, by going back to the simple 
case of iV = 1 and making the following plausible ansatz for the tachyon couplings: 

oo 

AS = ^2 fa e h( * +ix) + h e^-^ (5.13) 

k=l 

Next, assume that x is determined consistently by making the above term stationary: 

iL as=o (5 - i4) 

which leads to: 

t k e k<t >e ikx = t k e k *e~ ikx (5.15) 

or 

tk 

Now symmetry of the theory that interchanges t± with t\. The unique 

solution of the above equation that respects this symmetry is: 

tl = e ~ ix , tx = e ix (5.17) 

The first of these equations is a miniature version of the Kontsevich-Miwa transform, 
which is the best we can expect since we chose N — 1. 

Suppose we use the above solution for £i, and think of it as determining x. If we 
also leave t\ arbitrary (this procedure does not explain why we should do that) then 
we find that: 



AS ~ *fe eH e ~ lkX ( 5 - 18 ) 



k=l 



and t\ is determined by 



h = e~ ix (5.19) 

which is more or less the right story for the tachyon couplings of the Woo model at 
N — 1. A more complete analysis, particularly at general N, would hopefully shed 
some light on the non-dynamical nature of X in this model. 
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Interestingly there is another matrix model in the literature, the normal matrix 
model|38|], which is very similar to the Penner model but treats t, i in a manifestly 



symmetric way. It is defined in terms of a complex matrix Z that is constrained to 
be normal: 

[Z,Z^} = (5.20) 

and the matrix integral is: 

f , trf - v ZZ^ + 0- M)\ogZZ^ - v V?° ,( t k Z k + t k Z^)) 

Z(t,t) NMM = j[dZdZ^}e V ^ k=lK V 

(5.21) 

(A variant of this model exists for any finite radius R). The relationship of this 
model to the model deserves to be examined in detail, but we will not be able 
to do this here. Let us just observe that if we parametrise: 

Z = e*e lX (5.22) 

then the unperturbed part of the action depends only on <p and is similar to the 
unperturbed model. 

The normality constraint amounts to: 

[$,X]=0 (5.23) 

and the tachyon perturbation of this action is: 

AS = -v Y^ih e lk ^ +lX) + i k e lk ^- tX ^ (5.24) 

k=l 

We see that in this model, both the $ and X directions are represented by dynam- 
ical matrices, and the tachyon perturbation is similar to the one proposed above in 
Eq.(|5T3D. 

Indeed, the normal matrix model is supposed to be equivalent to the model 
at large N, as both are equivalent to c = 1, R = 1 by construction, though a precise 
equivalence in terms of matrices has not been demonstrated. In view of our conjec- 
ture that the action of the Woo model is that on N D-instantons with the transverse 
coordinate X treated as non-dynamical, it seems natural to speculate that the NMM 
describes the same system but with dynamical X. It is plausible that both mod- 
els become equivalent at large N, though superficially they do not appear to be 
equivalent at finite N. 

6. Conclusions 

We have seen that c=l,i?=lisa very special theory. It is completely solved to 
all orders in perturbation theory, and its free energy and correlators are known in a 
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genus expansion. Its partition function is related to the topology of moduli space, 
while the correlators are related to integrable hierarchies. 

All the above information is elegantly encoded in the matrix model. Here 
we have shown that this model maps to a matrix model with a Liouville-plus-linear 
interaction. This is proposed to be the matrix theory on iV D-instantons of the c = 1 
string theory. It is important to confirm whether this proposal is correct and/or 
needs modification. 

Extensions of these ideas and relationships from c = 1 to c = 1 (type OA, OB) 
strings will be interesting to pursue. In the supersymmetric case, there is naively 
no self-dual radius since the T-dual of type OB theory is type OA theory. However, 
there is a subtle "affine quotient" theory discovered in Refs.|64|, and studied in 
Ref. 0, which is T-dual to itself and therefore does have a self-dual radius. It is this 
theory to which the considerations of this article are most likely to generalise. 

Traditionally, topological string theories are tied to the genus expansion by their 
very definition, in that the partition function and amplitudes are defined genus by 
genus. But if, as is currently believed, the type OA and OB string theories truly 
have a nonperturbative completion, then the corresponding topological formulations 
might have one as well. If the latter can be summarised as matrix models, it is 
reasonable to hope that these super-analogues of the Penner and matrix models 
are not just generating functions of a genus expansion, but also contain physically 
meaningful nonperturbative terms. This could prove extremely important for the 
understanding of nonperturbative string theory. 

The above observations could also have a bearing on the many topological the- 
ories related to c = 1, R = 1, such as the topological black hole and topological 
conifold. 
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